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The present analysis is based on the complete first law flow equation i.e. it includes the kinetic energy term 
that was neglected in earlier papers. It is shown that the results obtained with the proposed equations, 
derived from an Entropy Maximizing Principle, agree quite well with reliable data published in the 
scientific/technical literature. 

The validity of the "theorem" of minimum entropy production (per unit time) was assessed in Refs 1 and 2 for 
the well documented case of single - phase fluid tube flow bifurcation from laminar to turbulent regimes. 
The analysis led to the following two main conclusions: 1) This theorem is not generally valid. 2) The stable 
type of flow regime is associated with a maximum specific (per unit mass) entropy change. More general broader 
statements of optimality in thermodynamics are now available in 3 . 

In order to avoid the range of Re in the so called "transition region" the analysis in 1, 2 was limited to Reynolds 
numbers one order of magnitude above the estimated critical Reynolds number R cr for the transition to turbulent 
flow and the kinetic energy term of the first law flow equation was neglected in line with normal engineering 
practice. 

In the present report the analysis is carried out down to Re„< Re cr and based on the complete flow equation i.e. 
it includes the kinetic term. As a result a new independent variable: the relative tube length or its inverse LID or 
D/L as well as two parameters: the laminar and turbulent kinetic energy factors M; and M t respectively, appear in 
the analysis and its solution. As the kinetic energy factors, or multipliers M play a key role in our analysis, a few 
words may be useful to clarify the definition of these parameters and provide a brief explanation on the methods 
adopted to assess their values under fully developed flow conditions for Newtonian fluids flowing in straight 
horizontal circular pipes. 

The kinetic energy term V/2g is only valid in the idealized world of inviscid fluid mechanics; due to the 
absence of friction the average velocity V also represents the local velocity that remains constant within the 
whole tube flow area. With real Newtonian fluids the local velocity is no longer constant and diminishes to 
zero at the pipe innerwall for both flow regimes. In the specific case of laminar regimes the actual local 
velocities could be determined accurately on the basis of Newton's laws of mechanics and from this a 
constant factor M ; = 2 determined to compute the correct kinetic energy term according to 2 (V 2 /2g). 
For turbulent regimes no exact theories are available to determine M, (as well / ( ) and we depend entirely 
on experiments. Fortunately M t ~ 1 and its impact on the prediction of Re cr is negligible compared to the 
laminar kinetic energy factor M; 

The same end results are obtained in 4 and the present report; the basic difference between these two commu- 
nications lies in the method adopted to predict the critical Reynolds number Re cr In 4 a quasi "empirical" 
approach was selected based on numerous plots produced with EXCEL. In the present report a purely analytical 
method was adopted based on a single premise: the extension of the validity of the Entropy Maximizing Principle 
proposed in 1, 2 to any Reynolds number, including the singularity at Re^ where D/L = O and LID — >oo, 

Although this topic may be of interest to some physicists, the theoretical prediction of the critical 
Reynolds number should prove particularly useful to scientists and engineers engaged in fluid mechanics 
and heat transfer mainly for several reasons that can be highlighted as follows. Soon after he demonstrated 
the existence of a "sinuous flow" regime and proposed an excellent estimate of Re cr = 2200 for smooth 
circular tubes, Reynolds derived what became known as the'Reynolds analogy", This shows that the internal 
hear transfer coefficient h for a circular pipe is directly proportional to the friction factor f t . Since the h t is 
usually much larger than hi it is advantageous to operate in turbulent regimes. Therefore an advanced 
knowledge of Re cr for instance in such radically new applications as micro-channel heat exchangers is quite 
important 4 . 
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Results 

Analysis. The well known 5 first law flow equation can be expressed as 
follows 

H = (L/D){V 2 /2g)\f+M(D/L)} (1) 

where H is the total head, V the average velocity and / the Darcy- 
Weissbach friction factor. Equation (1) is valid for laminar with//, 
Mf and / t , M t for turbulent flow. It is convenient to define a 
dimensionless relative head 

Hrd=H/Hfrict,<x, (2) 

where H frio,™ is the head loss due only to friction at Reco where the 
laminar and turbulent friction factor curves intersect or 

H frict , m = (L/D)(V 2 /2g)\f*_] (3) 

with 

U=f,=fiatRe=Re x (4) 

Since according to 1,2,6a , H = T As and Hf rictm = T Asy- nd)0O with T and 
As the absolute temperature and specific entropy change, it is clear 
that under constant temperature conditions 

H rd = AS rd (5) 

and in view of (1) to (5) 

H rd = AS re , = [f +M(D/L)]// C0 (6) 

The relative dimensionless entropy change As rd is used from here on 
instead of H re \- and it follows from (6) 

AS re lJ = \fl+Ml(D/L)]/f ao (7) 

AS reU = \f t +M t {D/L)]/f w (8) 

By treating both As re i <t and As re [j as independent mathematical 
functions, it is permissible to set 

AS rd J = AS rd ,t (9) 

yielding with (7) and (8) 

f t +M t (D/L)=f,+Mi(D/L) (10) 

noting that (10) is now independent of/^. and can be reduced to 

D/L=(f t -fi)/(M,-M t ) (11) 

It is interesting to note that as long as M; > M t and ( M; — M, ) is a 
finite number, then (11) yields D/L = 0, therefore LID — > 00 at Re m , as 
it should. 

If Mi and M t were exactly equal (11) would yield D/L = 0/0; 
however the two kinetic terms in (10) would cancel each other and 
the result would be:/ = f t = = f cr and _Re„ = Re cr . It is due to the 
fact that Mi I M t = 2 > 1 in the case of smooth circular tubes with 
fully developed flow at inlet that Re cr /Re« = 2.55 > 1 as shown in 4 . 

Minimum Re cn j d with fully developed flow at inlet of circular 
pipes. According to 5 M; = 2 and M t ~ I; introducing these values 
in (1 1) then the expression of D/L thus obtained in both sides of (10) 
yields the same equation for As re ;, ; and As re /, t 

AS rd = [2f t -f,} (12) 

As this mathematical function of Re has a maximum, we postulate 
that Re„ occurs there in order to maximize the entropy change 
according to 1 ' 2 . Therefore Re cr can be determined from 

dAS rei /dRe = 0 (13) 

or, with (12) 

d/,/dRe= (l/2)d/;/dRe= (M t /M ( )d/,/dRe (14) 



which means that the bifurcation to turbulent flow occurs at 
Re = Re cr where, in a linear coordinate system, the slope of f t is 
one half that of/;. 

Maximum Re cri f ud with fully undeveloped flow at inlet of circular 
pipes. This occurs when M ; ^- 1, for instance 1.01 or 1.001.... and 
M, = 1. Following the same procedure as in the previous case yields 
as a limit 

dft/dRe={l)df,/dRe={M,/Mi)dfi/dRe (15) 

The slope of / is identical to that of/ in this limiting case, meaning 
that Re„ t j ud must be significantly larger than Re cn f d . This is corro- 
borated to some extent by the following excerpts from the abstract of 
presented at the 2007 Annual Review of Fluid Mechanics: "Experi- 
ments on pipe flow,.... show that triggering turbulence depends 
sensitively on initial conditions." 

General case Re cn f d < Re cr < Re cr> ^ with Partially Developed 
Flow at pipe Inlet. It is easy to show that (14) and (15) are individual 
cases of the following more general equations 

AS rd = {Mif t -M t f,)/{M,-M t ) (16) 

d/ t /dRe=(M t /M,)d//dRe (17) 

for circular pipes: 

M, = 2^ = l,M f = l and M t /Mi= l/2-> = l 

It is noteworthy that the above equations were derived without 
using any explicit equations or correlations for/ or/, on the basis of 
an "Entropy Maximizing Principle", equation (13), and should there- 
fore be applicable with any reliable empirical correlation of / and 
valid for any pipe size, including micro tubes and channels. 

Prediction ofRe CTi f d with fully developed flow at inlet of smooth 
circular pipes. To validate the proposed theoretical equations it is 
necessary to compute Re cr with equation (14) and demonstrate that 
the critical Reynolds numbers obtained in this manner substantially 
agree with the most reliable data published in the scientific/technical 
literature. As the exact equation of/ in terms of Re is well known the 
right hand side of (14) presents no problem. For the left hand side we 
are faced with the following challenge: which one among the many 
empirical correlations proposed for the turbulent friction factor 
is best suited to our present purpose? Koo's excellent empirical 
correlation 7 based on his 1932 thesis at MIT is selected here for the 
following reasons. 

The values of/ obtained with four different correlations are shown 
in [5b]. Two of them, are clearly more accurate than the others, 
namely Petukhov's 9 and Koo's 7 correlations with the following 
Re validity limits 10,000 to 3,000,000 and 3,000 to 300,000 re- 
spectively. Since we focus in this study on Re below 10,000, Koo's 
correlation is the better choice; furthermore the maximum deviation 
between these two correlations being less than 1.5%, Petukhov's 
correlation is contained within the ±5% accuracy band of Koo's 
correlation based on 1380 most reliable test resuls. 

By introducing in (14) Koo's correlation for / and the exact 
equation of / for the laminar regime, we obtain d(0.0056 + o.5 
Re- °- 32 ) / dRe = 0.5 d( 64 Re~ 1 ) / dRe, yielding Re 0 63 = 200 resulting 
in Re cr = 2400. In order to make allowance for its ± 5 % accuracy, 
it is only necessary to multiply Koo's correlation by 0.95 and 1.05, 
resulting in a possible range of Re cr = 2160 to 2610. 

The uppermost value Re cr = 2610 < 3000, as it should, since this 
is the lower validity limit of Koo's correlation for turbulent flow 
based on tests with fully developed flow at tube inlet extending down 
to Re = 3000. 

The lowest value Re cr = 2160 ~ 2200, below which the flow 
remains laminar regardless of inlet conditions according to 5,6b . It is 
noteworthy that Osborne Reynolds, following his discovery of 
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Figure 1 | FANNO Flow Curve for Incompressible Fluids. 

Sinuous Flows'' in 1884, proposed a value of 2200 for what was later 
called the Reynolds Number! 

These positive results validate the proposed basic equation (17) 
with which they were obtained. 

Conversely, if one accepts this conclusion and remembers that 
(17) was derived from the application of the Principle of Entropy 
Maximization (13), then the principle itself has been validated. 

Estimate of Re cn ^ with fully undeveloped flow at inlet of smooth 
circular pipe. The impact of the laminar factor M ; on Re„ can be 
roughly estimated (as Koo's correlation based on tests with fully 
developed inlet flow is no longer as accurate) with (17) and 0.95./ t 
thus obtaining Re 06 " = 420 and Re cr ~ 7300. 

This estimate is sufficiently accurate to demonstrate the significant 
impact of M; on Re cr and to explain why Re cr > 3000 have been 
reported with smooth tubes as indicated in the 2009 Moody Chart 
shown in 5 . 

Discussion 

A most succinct definition 11 of engineering thermodynamics is: 
The Science of Energy and Entropy. Accordingly our analysis moved 
from conventional Fluid Mechanics to Classical Thermodynamics from 
equation (7) onwards and led, based on an entropy maximization 
principle, to the main equations (9), (11), (12), (13), (16) and (17). 

At some point it became clear to the authors that equation (16) 
would lead to an incompressible fluid F anno flow curve. Indeed all 
individual curves of As re ; in terms of Re produced with (16) can be 
reduced to a single common curve by changing the coordinates as 
follows: Replace As re i by As re ;/( As re i) max with (As re i) max obtained with 
(16) at_Re cr and substitute V/V cr = RelRe cr for Re. The curve obtained 
in this manner shown in Fig. 1 (reproduced from 4 ) is a FANNO 
FLOW curve according to the Extended Definition of Fanno Flow 
in 12 . The notes in green in figure 1 show the similarity with 
the more complex transition from sub-sonic to super-sonic flow 
regimes. 



One of the important features of classical thermodynamics is that 
its laws are independent of the size and geometry of the systems 
under consideration. This fact permits extending the application of 
(17) beyond smooth tubes to include non-smooth pipes and constant 
cross section channels of any shape and size i. e. including micro 
tubes and channels used in micro heat exchangers 3 , provided suitably 
accurate f and/; correlations as well as inlet parameters M; and M t 
are available with fully developed flow at pipe inlet. 

The need and means of suppressing the laminar flow regimes in 
order to reduce the size of compact tubular heat exchangers is illu- 
strated with the help of an internally finned tube design (INF) of 
rather complex geometry compared to standard circular pipes in 4 . 
This INF design patented in the 1940's was basically a forerunner of 
modern micro-channels heat exchangers. 

For scientists/engineers active in thermal sciences the following 
list of key words may be helpful to identify the novel or little known 
features included in this short report: Fluid flow thermodynamics; 
Maximizing entropy; Incompressible FANNO flow; laminar kinetic 
energy factor and; Impact of relative tube length. 
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